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Abstract 

Classical theory of nucleation based on Becker-Doering equations and 
coarsening for a binary alloy. 



1 Introduction 

The purpose of this chapter is to explain the classical kinetic theory of nucle- 
ation in a context simpler than polymer crystallization. Many theories start by 
assuming that polymer crystallization is an activated process involving crossing 
of a free energy barrier pp. The latter separates two accessible stable states of 
the system such as monomer solution and crystal. This general setting for acti- 
vated processes can be used to describe the formation of a crystal from a liquid 
cooled below its freezing point [2J , precipitation and coarsening of binary alloys 
colloidal crystallization chemical reactions [5], polymer crystallization 
El 0, etc. In all these cases, the theory of homogeneous isothermal nucle- 
ation provides a framework to study the processes of formation of nucleii from 
density fluctuations, and their growth until different nucleii impinge upon each 
other. In the early stages of these processes, nucleii of solid phase are formed 
and grow by incorporating particles from the surrounding liquid phase. There 
is a critical value for the radius of a nucleus that depends on a chemical drive 
potential, which is proportional to the supersaturation for small values thereof. 
In this limit, the critical radius is inversely proportional to the supersaturation. 
At the beginning of the nucleation process, nucleii have small critical radius and 
new clusters are being created at a non-negligible rate. As the size of existing 
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clusters increases, there are less particles in the liquid phase, the supersatu- 
ration decreases and the critical radius increases. Then it is harder for new 
clusters to spontaneously appear from density fluctuations. What happens is 
that supercritical clusters (whose radii are larger than the critical one) keep 
growing at the expense of subcritical clusters, that in turn keep losing particles. 
The size of the nucleii is still small compared to the average distance between 
them, so that impingement processes (in which two or more clusters touch and 
interaction between them dominates their growth) can be ignored. This stage 
of free deterministic growth is called coarsening UJ. 

A convenient framework to describe nucleation and coarsening is the classi- 
cal Becker-Doring kinetic theory. We assume that the dominant processes for 
nucleus growth or shrinking are addition or subtraction of one particle. Nu- 
cleation is thus treated as a chain reaction whereby nucleii of n particles are 
created by adding one particle to a nucleus of n — 1 particles, or subtracting one 
particle from a nucleus with n + 1 particles. We can then write rate equations 
for the number density of nucleii of n particles by using the law of mass action. 
The kinetic rate constants for the processes of addition and depletion have to 
be determined by using specific information from the physical process we are 
trying to model. Typically we impose detailed balance which implies that the 
ratio of rate constants is proportional to the exponential of the free energy cost 
of adding one particle to a nucleus of n particles (in units of fc^T). This leaves 
one undetermined rate constant. There are different ways of finding the missing 
constant. One way is to postulate a microscopic theory for particle interaction 
and use Statistical Mechanics to determine the free energy of a cluster [Hj. A 
different point of view is to impose that our rate constants should provide a 
description of coarsening compatible with the macroscopic description in terms 
of balance equations. We shall illustrate this second point of view and be led 
to a Smoluchowski equation from which the Lifshitz-Slyozov coarsening theory 
follows |10| . 

The structure of this paper is as follows. We present the Becker-Doring ki- 
netic equations for cluster with n particles in Section [21 One relation between 
the two rate constants of this theory follows from detailed balance. The other 
rate constant has to be determined by comparison with the known macroscopic 
equation for the growth of cluster radii. In the small supersaturation limit, the 
Becker-Doring equations can be approximated by a Smoluchowski equation for 
the distribution function of cluster radii. Its drift term yields the growth of clus- 
ter radii in terms of the missing rate constant. To compare with experimental 
data, we consider the case of coarsening of a binary alloy 3 . In Sections to 
we review phase equilibria, macroscopic kinetics of precipitate and matrix atoms 
and the quasistatic limit of the kinetic equations, respectively. As a result, we 
find the growth of the radius of a supercritical cluster in terms of macroscopic 
parameters. Comparison with the results in Section [21 yields the sought rate 
constant; see Section [S] Numerical values for all the parameters involved in our 
theories can be calculated from experimental data as explained in Section [SI A 
discussion of our results constitutes the last Section. 
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2 Kinetics of clusters 



Let us assume that we have two stable phases characterized by different values, 
c\ and C2 of the number density c. Phase 1 is solution and Phase 2 precipitate. 
Or Phase 1 is the liquid and Phase 2 the crystal phase. Initially all precipitate 
particles are in Phase 1. Classic Becker-Doring (BD) kinetics treats nucleation 
as a chain reaction whereby nucleii (assumed to be spherical) of n precipitate 
particles are created by adding one particle to a nucleus with n — 1 particles, 
or subtracting one particle from a nucleus with n + 1 particles. This chain 
reaction scheme is natural for the situation that BD has in mind, in which bulk 
precipitate phase consists only of precipitate particles so distinction between 
particles in nucleus or in solution is clear. 

Let p n be the number density of nucleii of n particles. The monomer density 
pi represents the concentration of precipitate particles in solution and as such it 
will be identified with the concentration of the macroscopic theory in Section 
|SJ Consider the reaction 

n + 1 ^ (n + 1). 

The forward reaction proceeds at a rate proportional to p\ p n with some rate 
constant k a . The backward reaction proceeds at a rate proportional to p n +i 
with rate constant kd- Hence the net rate of creation of (n + l)-clusters from 
n-clusters per unit volume is the flux 

jn = ka.nPlPn — kd.n+lPn+l- (1) 

The fact that the rate constants depend on cluster size has been explicitly 
indicated in JQ). Net rate of creation of n-clusters is due to their creation from 
(n — l)-clusters minus the rate of creation of (n + l)-clusters from n-clusters, 

Pn = jn-1 - jn = ~D- jn, n > 2. (2) 

This formulation specifies the evolutions of pi, P3, ■ ■ ■ with p± = Coo given. The 
number density of precipitate particles in n-clusters is np n and the density c 
(equal to the initial concentration of precipitate particles in the solution) of all 
precipitate particles is 

oo 

c=^2np n , (3) 

n=l 

or equivalently, 

oo 

c- ci = + 22 np n , (4) 

n=2 

where we have defined 7 = c — ci , 700 = Coc — ci . There is conservation of 
all precipitate particles so c or 70 = c — c\ are constant. (@J establishes that 
the constant initial concentration disturbance 70 is sum of the disturbance of 
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precipitate particles, 7oo, and the number density of particles in all cluster sizes 
n > 2. Given the constraint 10, the evolution of all p„, including pi, is specified. 

A most essential point of BD kinetics is identification of rate constants k a 
and kd- Ideally, this would be based on basic energetics and dynamics at the 
microscopic level, but a complete realization of this ideal is clearly elusive. Here 
is what is done: The ratio is determined by detailed balance. Equilibrium, if 
achievable, is described by a zero flux. Setting j n = in Q implies 

ka,n Pn+l /-\ 

(5) 



kd,n+l PlPn 

Standard equilibrium physicochemical theory states that 

Pn+l -Oil 



e 



(6) 



where p n is the free energy cost of creating an (n+l)-particle nucleus from an n- 
particle nucleus relative to the state of no nucleus, r = fc^T is the temperature 
measured in units of energy. Clearly, p n — G n+ i — G„ (G„ is the free energy of 
a n-cluster), so (0 becomes 



Pn+l 



e 



and JSJ now reads 



Q„ +1 -G» 

plk a ,n = e kd,n+l- 



Thus formula for the flux is now 

jn = kd, n +l |e T p n - pn+l t- (7) 

The equilibrium considered in the detailed balance argument is achievable 
only if G n — > +oo as n — > oo, whereas in a supersaturated solution, G n achieves 
a maximum for finite n and then G n — > — oo as n — > oo. The determination of 
the ratio k a /kd is assumed to hold regardless. 

What is known about G n l Microscopic models for G n (or, equivalently, the 
cluster partition function Q n = ^2 K e ub ^ K ^ T — e _G '*/ r , where U is the binding 
energy per pair of particles in the cluster of n particles, b{K) is the number 
of nearest-neighbor pairs of particles in the cluster K, and the sum is over 
all translationally inequivalent n-particle clusters) are described in [%1 1111 H2*] . 
We would like to follow here a simpler approach, consisting of identifying the 
resulting expressions for large spherical nucleii with known facts about radius 
growth in the quasistatic approximation. For nucleii of macroscopic size n» 1, 
n can be written in terms of the cluster radius a by 

n = —c 2 a i 1 (8) 
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and G n ~ G(a), where G(a) is the free energy of a nucleus of radius a as deter- 
mined by continuum theory. For nucleii of only a few particles, this asymptotic 
correspondence with continuum theory breaks down. But if the critical nucleus 
has n> 1 particles, the continuum approximation works for n on the order of 
critical cluster size. 

In the limit \(G n +i — G n )/r\ <C 1, formula J7J) for the flux reduces to 



where D±h n = ± (h n ±i — h n ). The basic evolution equation now reads 



This equation looks like a spatially discretized Smoluchowski equation. Asymp- 
totic replacement of difference operators D + , D- by derivatives is justified if the 
relative changes in G n , p n and kg when n increases by one are small. Here we 
will follow the simple procedure of formulating the continuum limit an checking 
its validity a posteriori. 

The space-like variable in I|1U|) is n. Experimental data usually contain 
histograms showing the distribution of nucleii in the space of their radii a, so 
we adopt the radius a as a more convenient space-like variable. The dependent 
variable should be p = p(a,t), the distribution of nucleii in space of radius a. 
Thus p(a,t) da is the number of nucleii per unit volume with radii in (a, a + da). 
Conversions n — > a, p n — > p are now determined. From |(HJ), it follows that the 
change da in a when n increases by 1 is given by 



In the general continuum theory of nucleii, the concentration of precipitate 
particles inside a nucleus, C2, is a function of the radius a. But for many 
experiments, deviations of C2 from its equilibrium value for a planar interface 
are negligible, so in l|llf) any term arising from a-dependence of C2 is dropped. 
p n is related to p(a, t) by 






(10) 



1 ~ 4-7T C2<2 2 da. 



(11) 




(12) 



Given any sequence h n with continuum approximation h{a), 




(13) 



It follows from (|12|l and l|13(l that the continuum limit of (|10|l is 




(14) 
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(Here kd is a function of a, to be specified). The constraint Q can be written 
as 



7oo+c 2 / — a 3 p(a,t)da = c - c\ = 70. (15) 
Jo 

As time elapses, it will be seen that the diffusive term in Eq. (|14fl becomes 
negligible in comparison with the drift term. The latter yields the following 
equation for radius growth: 

*" 9G (16) 



(47TC2a 2 ) 2 r 9a 

We now present a macroscopic theory that gives an explicit expression for d 
that can be compared to experimental data. Then kd{a) can be determined and 
this will specify the limit of kdn for large n. 



3 Phase equilibria of a binary material 
3.1 Phase equilibria of a binary material 

Let us consider a medium consisting of two different particles. The more abun- 
dant type is called "matrix" , the other "precipitate" . Suppose that we have 
a uniform mixture at fixed temperature and pressure. Let \x be the chemical 
potential , i.e., the free energy cost of adding one precipitate particle to a pre- 
existing solution. It is a function of the number density of precipitate particles, 
c: 

/i = /i(c). (17) 

Let us now derive the relationship of the chemical potential, /i(c), to the bulk 
free energy density, g{c). We shall add one precipitate particle to a solution of 
total volume V. Then the free energy changes from g{c)V to g(c)V + /i(c), but 
the volume changes from V to V = V + u(c), where v(c) is the specific volume 
of a precipitate particle in a solution of number density c. The number density 
changes from c to d = (cV + l)/(V + v), and therefore the new free energy is 
also expressed as g{d)V . Hence we get the identity 

g(c)V + rtc)=gOy^) (Y + u). (18) 

Since cV > 1, 1/ C V, this identity reduces to 

H(c)=g'(c) + v(c){g(c)-cg'(c)}. (19) 

One can just as easily consider the chemical potential /l(c), which is the free 
energy cost of adding one matrix particle to the solution. Adding one matrix 
particle changes the free energy to g(c)V + /Z(c). The volume of the solution 
changes now to V = V + v(c) , where 77(c) is the specific volume of a matrix 
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particle in solution. The concentration c of the precipitate changes to d = 
cV/(V + v), and therefore the new free energy is g(d)V. Hence we get an 
identity analogous to (jT%|) . 

g(c)V + /1(c) = g (t^—) (V + V). (20) 



y+v 

Again the conditions cV S> 1, V -C V lead to an asymptotic reduction of (|2H . 

/7(c) = 17(c) { ff ( c )-c S '(c)}. (21) 

3.2 Phase equilibria in a dilute solution 

Suppose that there are two stable phases characterized by different values, C\ 
and C2 of the number density c, and that these phases occupy adjacent half 
spaces separated by a planar interface. Imagine that one precipitate particle 
is removed from phase 1 and dropped in phase 2. The free energy cost is 
M c 2) — A*(ci). In equilibrium, the energy cost is to be zero, therefore 

M = A*(Pa) - A«(ci) = 0. (22) 

Similarly, the energy cost of removing a matrix particle from phase 1 and drop- 
ping it in phase 2 is also to be zero, so 

[p] = 0. (23) 

The possible existence of multiple phases is determined by the structure of g(c). 
Consider a dilute solution with c much smaller than the total atomic density. In 
this case, the specific volumes v and V of precipitate and matrix particles should 
be nearly constants independent of c: "crowding" effects should be insignificant. 
In this case, (|21|l and (|23|l imply 

[g(c)-cg'(c)} = [g-cg , ]=0. (24) 

Given this result, it now follows from (|T§|) and $Z2& that 

y\ = o. 

Hence, 

g'(ci) = g'(c2) — M (common value), (25) 

and 

[g] - [c] M = 0. (26) 

Given g(c) with g"(c) < in some interval of c, and g"(c) > outside, one 
discerns the standard geometrical construction of the solution to (|25|l and l|26|l 
for ci and c-i. This is depicted in Figure ^ 
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Figure 1: Geometrical construction of stable phase equilibria. 



3.3 Critical nucleus 

Consider a spherical nucleus of phase 2 surrounded by phase 1 . The energy cost 
of adding one precipitate particle to this nucleus is 

\p] + 8irar dr. 

Here r is the initial radius of the nucleus, dr is the change in radius due to 
adding one precipitate particle, and a is the surface tension. One has 

Ai:r 2 dr = v 



so the energy cost can be written as 

lav 



r 



For a nucleus in equilibrium, this energy cost is zero, therefore 

[,] = -—■ (27) 
r 

Now we add one matrix particle to the nucleus. No energy cost for this process 
implies 

[/!]=—-. (28) 
After substituting for ~p from l|21|) . this equation becomes 

\g-cg'] = -y. (29) 



Now substitute (JTSJl for (i and for 2a/r in (O to get 



[9'] = 0, (30) 

which is the same as in the case of a planar interface, l|25[) . Given the concen- 
tration ci of phase 1, this equation determines the concentration ci inside the 
nucleus, and then the Gibbs-Thomson relation (|27[1 determines the radius of the 
nucleus, r. This determination simplifies when the concentrations c\ and C2 are 
near "planar" values and 71 and 72 are deviations from the planar values. I|30[l 
together with 1|24[1 for the planar case imply 

[g"i] = o. (3i) 

Let us denote the common values of g"{ci) 71 and g"(c 2 ) 72 by to. The variation 
of [g - c g'\ in is 

-[cg" 7 ]=-[c]m. (32) 

Hence (gUl gives 

[c] to = — . (33) 
r 



4 Macroscopic kinetics 

4.1 Balance equations and jump conditions 

Let c(x,t), c(x,t) denote the macroscopic number densities of precipitate and 
matrix atoms. Local volume fractions of precipitate and matrix atoms are vc 
and vc, respectively. Since matrix and precipitate atoms fill space leaving no 
gaps, we have the space filling condition 

vc + Vc=l. (34) 

In conventional kinetics, the density of precipitate is locally conserved, with a 
flux proportional to the gradient of the precipitate chemical potential ^(c), 

Q + V • (— (5 V/i) = 0, 



or 

c t =V-(£)Vc), D = 8{c)n'{c). (35) 

Here 8(c) is a positive mobility coefficient and — SVfi = —DVc is the flux of c. 
This flux is formally a diffusion with diffusion coefficient D — 5(c) /i'(c). Given 
li(c) as in H19f) . and provided v and v do not depend on c (dilute solution), 

D = 5(c) fj,'(c) = 5(1 - vc) g"(c) = 5Vcg"(c). (36) 

In the last equality, the space filling condition has been used to replace 1 — vc by 
vc. For stable bulk phases, D must be positive, and |3*6^l then implies g"(c) > 0. 
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The description of matrix transport is essentially the same. The flux of matrix 
concentration c is 

-5Vp = -6jtfVc = -DVc (37) 

where ~p(c) is the chemical potential of the matrix atoms as a function of the 
matrix concentration c, 6(c) is the mobility of the matrix atoms, and D is the 
matrix diffusion coefficient given by 

D = 5ji'. 

Let g(c) be the free energy density as a function ofc. The space filling condition 
and 17(c) = 77(c) imply that 

M = 9'(c) +V{g(c) -cg'(c)}, 

which is totally symmetric to l|19|) . Then the diffusion coefficient D is related 
to 77(c) by a formula symmetric to (|36|l . 

D = 6vcg"(c). (38) 

The space filling condition leads to a relation between the mobilities 6 and 6. 
The linear combination vc + vc is locally conserved with flux —vDVc — vDVc. 
But vc + Vc = 1, so this flux is divergence free. Let C be any closed surface. 
Use of divergence theorem yields 



/ (v Dc n + v Dc n ) da = 0. 
Jc 



The space filling condition implies vc n — —vc ni so that we get \ c v(D — 
D) c n da = 0. This holds for all concentration fields c and closed surfaces C. 
Hence, 

D = D (39) 

and by lj35|> and (|35j> . 

6 (c) 77 c g" (c) = 5(c) vc g" (c) . (40) 

Now, 

5(c) = g(c) = g Q (1 - cv)^j , 
therefore, provided again that v and 77 do not depend on c, 

g"(c) = (=) g"(c), 

and (|40|l becomes 

(5iyc = <577c. (41) 
This is the relation between mobilities. 
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The integral form of Equation (|35[l informs the upcoming discussion of 
boundary conditions on a phase interface. Let R — R(t) be a time sequence of 
closed regions in which c = c(x,t) is a smooth solution of (|35|) . The number of 
precipitate particles in R is 



N = I cdx, 
Jr 

and the time rate of change is 

N = ctdx + Uc da, 
Jr JdR 

where U is the outward normal velocity of dR. Using i|35|) to substitute for c t 
above, and then the divergence theorem, we obtain 



N = [ (Uc + Dc n )da. 

JdR 



The interpretation of this equation is clear: The influx of precipitate atoms per 
unit area on dR is 

Uc + Dcn. (42) 

Suppose now that there is a region R\ of matrix surrounding a region R2 
of precipitate. The influx of precipitate atoms per unit area on interface C is 
given by l)42|) with c, c„ evaluated on the precipitate side of C. The outflux 
of precipitate atoms from the surrounding matrix into Ri is also given by l|42[l 
with c, c„ evaluated on the matrix side of C . Since precipitate atoms do not 
accumulate on C to form a surface density, the following jump condition holds 

[Uc + D c„] =0 on C. (43) 

In summary, conservation of precipitate is expressed by the diffusion equa- 
tion (|35|l and the associated jump condition i|43|) . The matrix density c satisfies 
the diffusion equation and jump condition with the same diffusion coefficient 
D. The space filling condition Ij34(l is automatically upheld by this kinetics. In 
addition there are "thermodynamic" jump conditions 

[</]=0, (44) 
[9 - eg'] = -2an, (45) 

expressing local equilibrium about the phase interface. These are in fact Equa- 
tions H29|) . (f 301) with 1/r replaced by the mean curvature k. Equations (|35() . 
143|l . (|44(l and (|45|l constitute a free boundary problem for the evolution of 
precipitate concentration c and the phase interfaces. 

4.2 Time evolution of Gibbs energy 

The total Gibbs free energy is 

G = Gi + G 2 + aS, (46) 
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where 



G x = [ gdx, G 2 = [ gdx (47) 

J Ri JRi 

are free energies of bulk matrix and precipitate phases, and S is the surface area 
of the phase interface, and a the surface tension. The time evolution of G under 
the kinetics of the free boundary problem (|35[) . 1)43 [I. I|44ll and (|45|) is examined 
here. To compute the rates of change of G\ and G2, it is useful to formulate a 
transport equation for the free energy density g(c): 

9t = 9' c t = g' V ■ (DVc) = V ■ (DVg) - D g" |Vc| 2 

or 

gt-V ■{DVg) = -Dg"\Vc\ 2 . (48) 

In stable bulk phases, g" > 0, therefore the source density in 148|) is generally 
negative. Let us now look at the rate of G2. 

G 2 = [ g t dx+ [ Ugda, (49) 
Jr 2 Jc 

where U is the normal velocity of the phase interface, positive if outward from 
precipitate. Inserting g t from l|48() and using the divergence theorem, (|49|l be- 
comes 

G 2 = [ (Ug + Dg n )\ 2 da - [ Dg"\\7c\ 2 dx. (50) 

JC JR 2 

In the surface integral, the subscript 2 means evaluation on precipitate side of 
interface. Similarly, 

G x =- [ (Ug + Dg^da- f Dg"\Vc\ 2 dx, (51) 
Jc JRx 

where subscript 1 means evaluation on matrix side of interface. The rate of 
change of the surface energy aS in 146(1 is given by the standard formula of 
differential geometry, 

aS = 2a [ nUda. (52) 
Jc 

Adding Equations l|50|l to (|52|) . we obtain the rate of change of the total free 
energy 

G= - [ Dg"\X7c\ 2 dx 

Jrx+R 2 

+ [ {[U g + D g n ] + IokU } da, (53) 
Jc 
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Here the jump [. . .] denotes values on precipitate side minus values on matrix 
side. Using the continuity condition [g'\ = 0, it follows that 

{Ug + Dg n ] = U\g}+g' [Dc n ]. 

In the right hand side, g' denotes a well defined value on the phase interface. 
By conservation jump condition (|43[1 . [D c n ] = — U [c], hence 

[Ug + Dg n ] = U[g-cg'\. (54) 

By the thermodynamic jump condition (|45|l . [g — eg'] = —2<jk, so finally, 

[Ug + Dg n ] = —2<7 f//c, 

and the energy rate formula (|53[) reduces to 

G = - I Dg"\X7c\ 2 dx. (55) 
JR1+R2 

The integral is negative definite. Notice that surface integral contributions to 
G cancel. This point is examined by direct physical argument to see what it 
really means. 

Recall that influx of precipitate atoms into precipitate phase per unit area 

is 

Uc + Dc n . 

The free energy of each precipitate atom changes by an amount [/i] = —2anv, 
according to 1|27[1 . as it crosses from matrix to precipitate. Hence there is a 
contribution to G of 

-2a(Uvc + Dvc n )n (56) 

per unit area of phase interface due to crossing of precipitate atoms. Similarly, 
influx of matrix atoms into precipitate phase per unit area is 

Uc + Dc„ 

and change in free energy for each matrix atom crossing into precipitate is 
[p] = —2ctkv, according to (|28|1 . Hence, crossing of matrix atoms gives another 
surface contribution to G, of 

- 2<j (UVc + DVc n ) k (57) 

per unit area. Adding l|5^jl and fBTjl yields surface contribution to G due to 
crossing of both types of atoms, 

— 2a {U (vc + 77c) + D (vc n + Fc„)} k = — 2aU n 

per unit area. Here the space-filling constraint has been used. From l|52|) it is 
seen that 2aU n can be identified as a rate of change of surface energy per unit 
area. Hence, the total rate of free energy production per unit area of phase 
interface is 

-2aUK + 2aUK = 0. 
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5 Quasistatic nuclei 



An isolated region R2 of precipitate phase, called a nucleus is assumed spherical, 
and concentration field c is assumed spherically symmetric. The kinetics is 
quasistatic if the time derivative in the diffusion equation 1(351) is negligible. Here 
kinetics is analyzed under the quasistatic assumption and regimes of validity are 
determined a posteriori by the criterion 

Ua aa „ . . 

-D = D « 1 < M > 

Here a is the radius of the nucleus, and the normal velocity U — a. a 2 / D 
represents the characteristic time of diffusive transport in the matrix phase 
surrounding the nucleus. The characteristic time associated with the kinetics of 
the radius is a/U — a /a. Kinetics is quasistatic if the time scale of the radius 
is much longer than the diffusion time in the surrounding matrix, as in (58} * 

Under assumptions of radial symmetry and quasistatic kinetics, the diffusion 
equation ((35(1 reduces to 

d r (r 2 Dc r ) = 0. (59) 
The conservation jump condition (|43(l reads 

a[c] = -[Dcr}. (60) 
The thermodynamic jump conditions (|44(l and 145(1 read 

w\ = 0, (6i) 

[9-cg'] = -^. (62) 

Given the value Coq of c as r — > 00, Equations 1(59(1 to ((62(1 determine an ordinary 
differential equation (ODE) for the nuclear radius a(t). The first integral of 1(59(1 
is 

r 2 Dc r = Q, or Dc r = Q. (63) 

Here Q is a function of time on r < a or r > a, but with a possible jump at 
r = a. Regularity of c at r = forces Q = on r < a. From ((63(1 it is now 
evident that 

[DCr]=-%, 
a 2 - 

where Q now refers to the value on r > a. The conservation jump condition 
l(6U|) now reads 



Q 

n2 



(64) 
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Figure 2: Geometrical construction of c(a— ) and c(a+). 



This is a differential equation for a(t), once the dependence of [c] and Q upon 
the radius are determined. The two thermodynamic jump conditions (|(jlf) and 
l|62l) determine c(a-) and c(a+), hence [c] as a function of a. Figure [21 shows 
the graphical construction of (|61fl and l|t)2|) . 

To determine Q, integration of (|63|l in r > a is required. Let ft.(c) be a 
function such that h'{c) — D(c). now reads 

d r h(c) = % , 
and integration from r = a to r — oo gives 

^(coo) - /i(c+) = — or Q = a{/i(coo) - Mc+)} . 
a 

Here c+ means c(a+). 164|) now reads 

a=fe)_M£±). (65 ) 
a[cj 

Since c+ and [c] are definite functions of a as determined by the thermodynamic 
jump conditions (|61() and 165|l is the required ODE for a(t). 



5.1 Small supersaturation 

A standard limit called small supersaturation is realized when the concentration 
is close to planar equilibrium values c = c\ in matrix phase r > a and c = C2 in 
precipitate phase r < a. Let 7 denote disturbance of c from planar equilibrium 
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values: 7 = c — c\ in matrix phase, 7 = c — C2 in precipitate phase. For 
|7| <C ci,C2, (|63j) reduces to 

a n . 66) 

a[c] 

Here D\ = D{c\) and [c] = C2 — Ci. Also 700 = j(r — 00) and 7+ = 7(0+). 7+ 
is determined from asymptotic limit of (|61|) and l|62|) : Reduction of l|61(l is 

Wl] — </'( c i)7+ = .9"( c 2)7- = (common value), (67) 
and given m, reduction of (|62|) is 

cm=— or j ' c 7+ = - . (68) 
a a 

Substituting this result for 7+ into 166|) gives the reduced ODE 

fll 1 / „ 2(7 \ 

'Sl7ao--n • ( 69 ) 



[c] gj' a \ " " a [c 



This equation indicates that clusters whose radii are smaller than the critical 
value 

a c = — tt , (70) 

9" loo [C\ 

shrink and disappear. Supercritical clusters of radius larger than the critical 
radius a = a c grow steadily according to (|7Tj|l . 

Is this small supersaturation kinetics consistent with the quasistatic criterion 
H58D ? Natural unit of a is a c given by l|7U|). which is the standard formula for 
critical radius in small supersaturation limit. Given a — 0(a c ), an order of 
magnitude estimate of a based on (|p^|) is 



a = O 

which can be rearranged as 



£i7o 



The analysis here is based on I7I <C Ci, C2. But for quasistatic kinetics, we need 

|7oo| « [c] =c 2 -C!. (72) 

The reduced ODE (|69|l indicates natural units of 7, space and time. The 
limit (|72|l is embodied by measuring 7 = c — c e in units of e c e , where e > is a 
gauge parameter and limit e — * 0+ is considered. Given this unit of 7, the order 
of magnitude of a c in (|70|l is l/e, I = cr/([c]<?"c e ). l/e is adopted as unit of length. 



16 



Finally, the unit of time which gives a as in lf7T|) is r e /e 3 , r e = (I 2 / ' D{) ([c]/c e ). 
This system of units is summarized in the following Table: 



Variable 


Unit 


7 EE C - C e 


ec e 


X 


l 1 — cr 
«' ~~ Wff"c c 


t 




G 


f 2 



Scaling Table for small super ■saturation. 



Given these units, the nondimensional version of (|69|l is 

a = - (700 J • (73) 



5.2 Quasistatic energetics 

Given evolution of precipitate concentration field c, corresponding changes in 
total free energy of medium are quantified by the rate formula (J5SJ. Now suppose 
concentration field c corresponds to a nucleus undergoing quasistatic evolution 
as set by criterion (|58|l . It seems reasonable to approximate the integral in l|55JI 
using the quasistatic approximation to c which satisfies Equations (|59|l to (|62|) . 
The rate formula l|55|) reduces to 

/>oo 

G = — Dg'clAnr 2 dr. 
Jq 

/From (|63|l . r 2 Dc r = Q for r > a, and r 2 Dc r — for r < a, so this reduces to 

/>oo 

G = -4^Q / g"c r dr = -4irQ { 5 '( Coo ) - g'(c+)}. 

Substituting for Q from l|64|) . 

G = -[c] { 5 '( Coo ) - .g'(c+)}(47ra 2 d). (74) 

In the right hand side, [c] and g'(c+) are definite functions of radius a, deter- 
mined by thermodynamic jump conditions (|61|l and (|62|l . Hence, in quasistatic 
evolution, G is effectively a function of the radius a. ^From (|74|l . it is seen that 
G(a) obeys the differential relation 

dG - -[c] {g\ Coo ) - g'(c + )}(4ira 2 da). (75) 

In conventional descriptions of nucleus energetics, the Gibbs free energy cost of 
a nucleus regarded as a function of the instantaneous radius a, independent of 
past history. The specific formula is 

Air 

G = Airaa 2 - g b — a 3 (76) 
o 
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or in differential form, 

dG = 8iro~ada — gbAira 2 da. (77) 

Here a is surface tension and 4iraa 2 represents surface energy. gt> is a constant 
with units of energy density sometimes called "chemical driving force" . It is 
the free energy released per unit volume of increase of precipitate phase. In its 
present form l|75|l does not have obvious correspondence to l)77p. A correspon- 
dence is brought out by reformulation of (|75|l with help of thermodynamic jump 
conditions l|HT|) and \<o2\. By IpxT]) . it follows that g'(c_) = g'(c+) = common 
value m, and hence gives 

[g] - [c]m= - — , 



or, equivalently, 



[c]g'(c + ) = [g] + - 



Hence (|75H becomes 

dG = 8iraada + {[g} - [c] g'(c 0O )}(47ra 2 da). (78) 

The first term on RHS is differential of surface energy, same as in (|77(l . An 
apparent correspondence between (|77|l and l|78(l is completed by identifying the 
chemical driving force gb, 

g b = -[g} + [c}g'( Coo ). (79) 

Recalling that [g] and [c] are functions of radius a as determined by thermo- 
dynamic jump conditions, it is evident that the chemical driving force l|79|l is 
generally a function of cluster radius, and not a constant as in conventional 
nucleus energetics. 

We can now write Eq. I|fi4[) for a in terms of G a = — {g'(coc) — g'{c+)} 47ra 2 [c] 
given by 1)75(1 . The result is 



f^D(c)dc Ga 



g'( Coo )-g'(c + ) 4^[c] 2 a3- 

5.3 Energetics at small supersaturation 

In the small supersaturation limit, J c c °° D(c) dc ~ D(c\) (c+ — Coo) — D\ (7 
Too), g'(coo) - g'{c+) ~ g'{ (7+ - 700) and (E9 reduces to 

. = G a D 1 
a 47r<[c] 2 a 3 ' 



(80) 



(81) 



Inserting the approximate chemical driving force g^ = [c]<7'(coo) — [g] ~ <7i7oo 
from l|79|l into the free energy (|76|l . we obtain 

G ~ iiraa 2 - [c] g'{ 7oo ( y . (82) 
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Notice that the asymptotic chemical driving force, 



9b ~ [c]g'{ Too, (83) 

is in fact a constant independent of radius a, as in the conventional wisdom. An 
alternative expression for gt is useful: In certain experiments surface tension a 
and critical radius a c are measured observables, so it is convenient to represent 
gi, in terms of a and a c , 

<* = -. (84) 

This formula follows from the condition G'{a c ) = 0. The natural unit of free 
energy is al 2 . This is entered into the last column of the scaling table. Then 
the dimensionless version of the free energy formula 182|) is 

47T 

G^Aira 2 - loo —a\ (85) 



5.4 Identification of rate constant in BD kinetics 

We can now compare Eq. for the growth of a (large) cluster radius in 
BD kinetics with the corresponding equations l|8U|l and 181|) obtained from our 
macroscopic description. We hnd 

t r°°D(c)dc / C9 \ 2 
S'(coo) ~5'(c+) \[c]J 

~f(§)%™, <sr) 

in the small supersaturation limit. 

With this determination of kd, the continuum limit equation (|14|l of clus- 
ter kinetics is completely specified. It is easy to check that the corresponding 
microscopic rate constant determined by Penrose et al ^2 IS] is a ls° propor- 
tional to a for large clusters. One might wonder about the micromolecular basis 
of 1)8 7|) . While that requires further work, here is a curious observation which 
might become relevant to this question: Recall the mobility 8 defined in the 
formulation of the macroscopic transport theory of Section 21 It is related to 
the diffusion D by ij^jl . 

D = S{l-vc) g". 
Hence the ratio D\ / g'{ in (|87|l is given by 

^=5 x {l-vc x ). (88) 
Si 

In the next Section, we show that the volume fraction vc\ of precipitate in 
matrix phase is small, vc\ ~ (0.04521 nm 3 ) (2.24 nm -3 ) w 0.10, for coarsening 
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experimental data in binary alloys Then Dx/g" » 8\ and formula (|87|l for 
fed reduces to 




Let us notice that this rate constant k^ is linear in the cluster radius so 
that it scales as ns with cluster size. This contrasts with the usual Turnbull- 
Fisher rate constant that scales as ns but it agrees with the microscopic 
considerations of Penrose et al m^]. The scaling 113 has been shown to yield 
the Lifshitz-Slyozov distribution function for cluster radii The latter is a 
roughly adequate description of coarsening 0E]- 



6 Material and energy parameters of kinetic the- 
ory determined from Xiao-Haasen data 

6.1 Small supersaturation 

The nucleation in Xiao-Haasen's (XH) paper takes place under low super- 
saturation: The initial sample has uniform composition, with mole fraction 
X = 0.12, or 12 % of Al in a Ni matrix. Equilibrium mole fraction of Al in 
matrix phase at annealing temperature of 773 K is xi = 0.101, or 10.1%. Equi- 
librium mole fraction of Al in precipitate phase is \i = 0.230, or 23 %. Hence, 
supersaturation as a function of equilibrium concentration has initial value 

X-Xi „, 0-120 -0.101 
Xi 0.101 

Mole fractions are converted into number densities: XH report a molar volume 
of precipitate phase V m « 27.16 x 10~ 6 m 3 . Conversion to an atomic volume 
by Avogadro's number gives 

Vm = ^i« 4 .51x 10- 29 m 3 

or v m « 0.0451 nm 3 . 

XH also report a lattice constant of a w 0.356 nm for the precipitate phase, 
and atomic volume corresponding to this lattice constant is v m — a 3 ss 0.0451 
nm 3 . It is clear that the molar volume V m was derived from the lattice con- 
stant. A lattice constant for the matrix phase is not reported explicitly, so it is 
presumably close to the value a w 0.356 nm of the precipitate phase. It seems 
there is an implicit assumption: Local structure of alloy in a lattice, with sites 
that can be occupied by Al or Ni atoms. In this case, atomic volumes of Al and 
Ni are de-facto the same, i.e., 

Vm = v m ~ 0.0451 nm 3 . 
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Now number densities of Al and Ni easily follow. For instance, ci, the equilib- 
rium number density of Al in matrix phase is 

Xi 0.101 , „„, o 

ci = — « nm -3 ps 2.24 nm -3 . 

0.0451 

Table 1 gives initial concentration of Al in matrix phase, and equilibrium con- 
centrations ci and C2 of Al in matrix and precipitate phases. 



c 


Cl 


C2 


2.66 


2.24 


5.10 



Table 1: Number densities (nm 3 ). 

6.2 Nucleation energetics 

In XH, the nucleus energy takes the classic form 

G = Altera 2 - g b ( y a 3 \ . (89) 

Here a is surface tension. A value a ps 0.014 J m~ 2 = 1.4 x 10~ 20 J nm -2 is 
deduced from interpreting coarsening data with the Lifshitz-Slyozov (LS) theory. 
XH deals with chemical driving force gb in two ways: 

(i) "Experimental" . The distribution of nuclei in the space of their radii goes 
through a transient phase with two peaks, separated by a local minimum at 
about 1.2 nm. XH conjecture that the initial radius a c is in fact this 1.2 nm. 
Estimate of gb now follows from (|83|l . 

g b ~ — ps 2.33 x 10~ 20 J nm" 3 . 
a c 

Given an experimental estimate of the critical radius, a c ps 1.2 nm at outset, 
one can estimate the number of atoms in the critical nucleus, both Al and 
entrained Ni: 

1 ^ al « 160. 

v 3 c 

Of these, 23% are Al, so there are 

n c = 0.23 x 160 = 37 

Al atoms in the critical nucleus. It seems that the critical nucleus is "just big 
enough" so energetics based on continuum theory applies. 

(ii) "Theory". In standard theories, gb is computed from both thermodynamic 
properties of precipitate and bulk phases. As such, it comes out as a constant 
independent of nucleus radius a. These derivations do not face up to the fine 
points of the real situation, summarized in the formula (|79|l for g b . So our 
approach is to stick with the determination of gb based on a and a Cl 

gb ~ ~ 2.33 x 10~ 20 J nm" 3 , 
a c 



21 



and then see what can be said about g(c). In l|83[l one can determine g'{ because 
all the other quantities are known. In fact, one gets 



9b 2.33 x 10" 20 J nm in _ 20 _ 3 

— ~ ~ 1 Q4 y 1(1 T nm 

c] 7oo (5.10 - 2.24) nm" 3 (2.66 - 2.24) nm" 3 



The annealing temperature of 773 K defines a basic unit of energy, 
r = (773 K) (1.38 x 1CT 23 J/K) » 1.07 x 10~ 20 J. 

One now has 

^«1.81nm 3 . 
r 

This is just one number imposed upon free energy function g(c) by the XH data, 
but it is sufficient to establish the 

Nonideal character of Al solution in matrix phase. 

Suppose the solution is ideal. Then the chemical potential of an Al particle 
in matrix phase is given by 

H(c) =Mi+Tln-, (90) 

Cl 

where /ii is the chemical potential when c = C\ is the planar solvability. Now 
the relation between /t(c) and g(c) is given by i|19|) . which is repeated here for 
easy reference, 

t i(c)=g'(c)+v(g-cg'). (91) 

^From this equation is evident that g"(c) gives information about //(c). In fact, 
differentiation of l|91|l yields 

//(c) = (1 - vc) g"{c) =*> /ii = (1 - i/ci) gf" . (92) 

Numerical value of ///r based upon previous value of g'{ turns out to be 

— w {1- (0.0451 nm 3 ) (2.24 nm" 3 )} (1.81 nm 3 ) « 1.63 nm 3 . 
r 



If the ideal solution formula l|90|l were correct, one would get 

0.45 nm 3 



t ci 2.24 nm -3 

which is 1 /274 of value that follows from XH parameters. That the solution of Al 
in Ni phase is not ideal was already known to XH. They in fact considered that 
our chemical driving force gb is sum of two terms: (i) a chemical driving force 
estimated from the activity of Al component at the concentrations \ and Xi> 
and (ii) and the elastic strain energy per unit volume. With the corresponding 
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expressions, they obtained a value for the critical radius, a c ss 1.7 nm, which is 
not too far from the experimental value, a c w 1.2 nm 

The experimentally derived values of surface tension a and chemical driving 
force gb in l|89|l set important parameters for macroscopic nucleation theory, 
namely: Energy barrier for nucleation, and typical free energy cost to add one 
Al particle to a nucleus. 

Energy barrier is given by 

Gnuc = G(a c ) = — aa\ w — (1.4 x 1CT 20 J nm" 2 ) (1.2 nm) 2 w 8.4 x 1CT 20 J. 
Energy barrier in units of thermal energy is 

^ „ 7. 8 . 

r 

A reasonable looking number. Notice that exponential 

e"^ « 3.7 x 10- 4 , 

which appears in the nucleation rate is not too small. This makes anthropo- 
morphic sense: In XH experiment, nucleation kinetics unfolds in hours and days 
time scales, not unduly taxing to humans. Evidently, the annealing temperature 
is tuned so as to achieve a "reasonable" nucleation rate. 

Free energy cost to add one particle. 

The number n of Al particles in nucleus is related to radius a by 

n= y c 2fl ^ a =(— J . (93) 

Substituting (|93(l for a in Ij89(l gives nucleus energy as a function of n, 

G n = (367r)3 ac 2 3 ni — g b c^ 1 n. (94) 

Free energy cost to add one particle to nucleus of n particles, in units of thermal 
energy r, is 

(J, n G n +i — G n . -i -2 i ,. .2 2 _j i 

— = — - = (367r) a ac 2 3 t {(n + l) 3 - n 3 } - gb c 2 r \ 

t r 

Since this formula is based on continuum theory, its validity requires n 3> 1, in 
which case it reduces to 

^ ~ l(36n)i ac'' 3 r^n-i ~ gtc^ 1 t- 1 . (95) 

T 3 

Substituting XH parameter values in RHS, 

^ ~ 1.43 n'i -0.42. (96) 
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In the limit n — > oo, we get |/i„/r| ~ 0.42. While less than 1, one would 
not call this value "small compared to one". For n — 37, corresponding to 
a critical nucleus, of course one gets /U n / T = 0. Hence there will be a range 
of n about n — n c = 37 in which \fi n / T \ C i In particular, |/in/ T l < 0.2 
in the rather generous interval 12 < n < 300. In this range of n, asymptotic 
reduction of discrete kinetic models such as Becker-Doring to a Smoluchowski 
partial differential equation (PDE) should be reasonable. 

7 Discussion 

In coarsening experiments, one starts from a situation of equilibrium at high 
temperature in which most clusters are monomers. Then the temperature is 
lowered to a value below the critical temperature, and kept there. Clusters 
are nucleated and grow, supersaturation changes with time so that nucleation 
of new clusters becomes unlikely, and the coarsening of clusters proceeds. As 
explained by Penrose et al m E]> this process is reasonably well described 
by the Becker-Doring model (better than by the Lifshitz-Slyozov distribution 
function), provided the volume fraction of precipitate is small. Let us describe 
the nucleation and coarsening processes in typical experiments such as XH and 
which parts thereof arc mathematically understood. 

The nucleation process described by the BD equations starts at t = with 
some initial value of p\ = and no supercritical clusters. According to the 
XH data, the energy barrier corresponding to the initial value of the critical 
nucleus is relatively high, Gnuc/ T ss 7.8, so that we may consider the clusters 
below critical size (n < n c ) to be in a quasistationary state. The flux across 
the energy barrier is then uniform and it supplies the source for coarsening of 
clusters larger than the critical size. As explained in Sectional there is a range 
of sizes (about the critical size) for which we may approximate the discrete BD 
kinetics by a continuum Smoluchowski equation for the distribution function p. 
The latter will describe the coarsening process and it should be approximately 
solved with a boundary condition obtained by matching to the solution of the 
BD equations for n < n c . For t > 0, supercritical clusters are created at the rate 
j per unit volume given in Eq. (|102|l below, and p\ starts to decrease. A small 
change of pi, 0(l/n c ) = 0(e 3 l 3 /c2), induces an O(l) relative change of j. There 
is a transient situation during which p becomes a bimodal distribution function 
with peaks at sub and supercritical sizes. As time evolves, the supercritical peak 
increases at the expense of the subcritical peak, which disappears given enough 
time. Then the resulting unimodal distribution evolves toward a function with 
the LS scaling. 

Currently it is known that the LS distribution function [H] is a solution of the 
Smoluchowski equation for a very special boundary condition at small cluster 
size PO]- Although the stability properties of the LS distribution function are 
not completely elucidated, it seems clear that the Smoluchowski equation may 
have other stable solutions that may match the quasistationary distribution 
at small cluster sizes. The appropriate solution of the Smoluchowski equation 
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should then describe the transient stage of coarsening. As the time advances, the 
peak of the distribution function at subcritical sizes decreases and disappears 
while the peak at supercritical sizes takes over. The latter should have the LS 
scaling to explain experimental |3] and numerical data [111 I12| . To carry out an 
asymptotic analysis of nucleation and coarsening providing the same qualitative 
description sketched here is a challenging future task. 
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9 Appendix: BD kinetics for n < n c 

The quasistationary state is a solution of the BD equations characterized by 
uniform flux, 

]n = k d ,n+l je - p n — Pn+l j = j, (97) 

for n < n c = C247raj?/3. At high temperature, before the experiment starts, we 
have the following equilibrium solution 

G n 

Ce r G n — Gi 

Peg-n = ~ce - . (98) 

To write the above approximation, we have assumed that G%/t G\/t 1 
and that G n increases with n. After nucleation and coarsening start, we shall 
assume that p n is close to its equilibrium value [given by the approximate expres- 
sion (|98|l at the correct temperature r], as n <C n c . Thus p n — 0(ce~( Gn ~ Gl ^' r ) 
if n < n c , and much smaller than this order if n > n c . This means that 
e - ^ Pn/c = 0(e Gl / T ) if n < n c , and that p n /c = o{e Gl / T ) if n 3> n c . 
Equation (|97|l can be written as 



g, 



+1 



e T Pn+i - e 



kd,n 



+ 1 



and therefore easily integrated under the condition e " p n — > as n — > oo: 

oo G i + i 

Pn=]J2 J— ■ (") 

The terms in this sum are largest for / ~ n c , at which Gi is maximum. For such 
integers, the continuum approximation holds, and we can write 

G i + 1 G "e OO 

£_^„f_^ e -^ (-«-.)». (100) 

' ' l—n 
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We have used G;+i — G„ c ~ —Ai:a(ai — a c ) 2 , for I + 1 close to n c . We now 
approximate ai = a c + X\J a jr in i|99|) . and 1 ~ 4TTC2a 2 dai = \-kcic? c \J t j a dx , 
so that lf§9")) becomes 



e ^ ^ ~ J_ M Vi 2 / e^'dx. (101) 



* (a-o c ) 



The equilibrium solution of the BD equations is l|98(l . If we impose that p n ~ 
p egi „ as n <C n c , (|9*S|) and (|101|l yield 

This constant flux is exponentially small because(G„ c — Gi)/r ~ G(a c )/r ^> 1. 
Notice that it is also proportional to the supersaturation 700. It is clear that a 
small change in the supersaturation, J7 = 0(7oo/"- c )j produces an O(l) change 
in n c and in G nc , Sn c = — 3n c ^7/7 00 and SG„ a — — 2G„ c <57/7 00 , and hence a 
significant relative change of j in i|102|) : 

« - expfffl!^^ • (103) 



J 

Notice that, in the continuum limit, the flux j n becomes 



g(°) 



D±t e t d 

jn ~ 



g'{[c] 2 Ana 8a V a 2 ) 

5 i'[c] 2 Ana V™ 2 da 8a \a 2 J J { ' 

The relation between drift and diffusion coefficients here is G a /r in agreement 
with the formulas provided by Nonequilibrium Thermodynamics; see Ref. UJ. 
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